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ABSTRACT 


The grand-canonical version of the thermal model is used to analyze the ratios of particle 
abundances measured in ultra-relativistic heavy-ion collisions. Exactly the same model is 
applied to study the heavy-ion reactions at BNL AGS, CERN SPS, and BNL RHIC. A very 
good description is achieved for Pb + Pb collisions at CERN SPS and for An + An collisions 
at BNL RHIC. In these two cases the value of the temperature characterizing the chemical 
freeze-out is practically the same: we find T chem = 168 ±3 MeV at SPS, and T chem = 165 ±7 
MeV at RHIC. On the other hand, the particle ratios measured in the collisions of lighter 
nuclei are described only in the qualitative way. We discuss also the effect of the possible 
in-medium modifications of hadron masses and widths on the thermal fits. For Pb + Pb 
collisions at CERN SPS and Au + Au collisions at BNL RHIC, we find that the y; 2 fits favor 
slightly a moderate, ~ 20%, decrease of the masses. In this case, the fits with the modified 
masses yield modified values of the optimal temperature and the baryon chemical potential. 
In-medium modifications of the widths have little effect on the fits, unless they are increased 
by a factor larger than 2. We study in detail the thermodynamic conditions characterizing 
the chemical freeze-out. In particular we find that the average baryon energy at freeze-out 
is 1.6 GeV, and the average meson energy is 0.9 GeV. This difference reflects a different 
behavior of the mass spectra of baryons and mesons. Similarities and differences between 
our calculations and other studies are thoroughly discussed. 
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Chapter 1 
Introduction 


1.1 Ultra-Relativistic Heavy-Ion Collisions 

The field of the ultra-relativistic heavy-ion collisions connects the nuclear physics with the 
elementary particle physics. In the high-energy particle physics the interactions are de¬ 
rived from first principles ( local gauge theories) and one deals with single particles ( leptons, 
hadrons, or quarks and gluons ). On the other hand, in the traditional nuclear physics the 
strong interaction is described by effective models, and the matter consists of extended 
complicated systems (nuclei). A unifying aspect of the high-energy nuclear collisions is an 
attempt to analyze the properties of dense hadronic matter in terms of elementary interac¬ 
tions. The fundamental theory of strong interactions, Quantum Chromodynamics (QCD), 
predicts that at high energy density, hadronic matter will turn into a plasma of deconfined 
quarks and gluons (QGP). The search for such a phase transition is the main motivation 
for the continuous experimental and theoretical efforts in the field [1-8]. One of the goals 
of the experimental program is to recreate (on a microscopic scale) the physical conditions 
that are thought to have existed in the early universe. Such astrophysical aspect of the 
ultra-relativistic heavy-ion collisions indicates once again on the interdisciplinary character 
of this new branch of physics. 

The experimental studies of the high-energy nuclear collisions started in 1986 at the 
Brookhaven National Laboratory (BNL) and at the European Organization for Nuclear 
Research (CERN). The Alternating Gradient Synchrotron (AGS) at BNL accelerated 28 Si 
beams (at 15 GeV per nucleon) and later 197 Au beams (at 11 GeV per nucleon). At CERN 
the Super Proton Synchrotron (SPS) delivered 16 0 and 32 S beams (at 60 and later at 200 
GeV per nucleon), which were followed in 1995 by 208 Pb beams (at 158 GeV per nucleon). 
The new era in the field started in 2000, when the Relativistic Heavy Ion Collider (RHIC) 
started operation at BNL. In the first run of RHIC, Au on Au reactions were studied at the 
center-of-mass energy yfs = 56 A GeV and yfs = 130 A GeV. In 2001, during the second 
run the full collision energy, yfs = 200 A GeV, and the full luminosity were achieved. 
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Over the last 15 years, the large amount of data has been accumulated. In particular, the 
data indicate that the hadronic matter at freeze-out (i.e., at the moment when the hadron 
interactions cease and the particles freely stream from the collision point to the detectors) is 
well described by the equilibrium distributions [9-23]. In addition, as one moves up from SPS 
to RHIC energy, approximately the same temperature and a significantly smaller baryonic 
chemical potential is observed in the central rapidity region. This fact constrains different 
possible models of the particle production. Actually, there is no single space-time model 
describing the whole collision process, since different degrees of freedom are important at 
different stages of a collision. As a consequence, each stage requires a different treatment: 
The initial moments are described with the help of partonic degrees of freedom (quarks and 
gluons). The intermediate stage is described typically in the framework of the relativistic 
hydrodynamics (in this case the degrees of freedom are collectively included in the equation 
of state). Finally, the last stage is described in purely hadronic models. 


place 

projectile 

target (s) 

E lab /A [GeV] 

Vlab 

y/sfA [GeV] 

BNL AGS 

Si 

Au, Pb 

14.6 

3.4 

5.4 

CERN SPS 

S 

Au, W, Pb 

200 

6.1 

19.4 

CERN SPS 

Pb 

Pb 

158 

5.8 

17.3 

BNL RHIC 

Au 

Au 

( 9000 ^ 


130 

^ 21300 J 

l 10.7 ) 

200 


Table 1.1: Heavy-ion reactions studied in this paper. 


1.2 Aim of this Work 


In this work we discuss the properties of matter at freeze-out, thus we are not interested 
in the microscopic mechanism of particle production, such as decaying color strings |24 
or a parton cascade f25|. We assume, however, that the real microscopic mechanism leads 


finally to creation of a locally equilibrated hadron gas. We note that in our approach one 
cannot draw any direct conclusions concerning formation of the quark-gluon plasma at the 















earlier stages of a collision. Such information is lost in the thermal equilibrium, which has 
no memory about earlier times. 

Our main purpose is to perform a thermal analysis of the particle yields. We are going to 
check whether the measured particle yields may be explained in a model which assumes full 
thermal and chemical equilibrium of the hadronic matter at freeze-out. In the last years such 
thermal approach has been used by several groups to study different types of reactions. The 
results of such investigations show that the thermal description of the particle production 
is quite successful. We have to have in mind, however, that different groups use different 
implementations of the model, and these implementations vary for each particular reaction. 
In this situation, the conclusions drawn from the vast applications of the thermal models 
may be not completely consistent. In this paper, to minimize such effects, one formulation 
of the model is used to describe most of the available data. In this way, exactly the same 
thermodynamic features of the particle production are used to characterize different colli¬ 
sions. This allows us to observe similarities and differences in the thermodynamic behavior 
of various colliding systems. 


The second and independent aim of this work is to include the possible in-medium mod¬ 
ifications of hadron properties into the thermal approach. In the scenario with the two 
different freeze-outs, i.e., with the chemical freeze-out preceding the thermal freeze-out, the 
study of the particle ratios reveals the information about the chemical freeze-out: the con¬ 
cept of the chemical freeze-out refers to the point at which inelastic collisions cease and all 
particle ratios are frozen, whereas the concept of the thermal freeze-out refers to the stage 
when all elastic collisions cease. If the hadronic matter at the chemical freeze-out is very hot 
and dense, we may take into account the mass and width modifications of hadrons. Such 
modifications are predicted by the effective theories of QCD [26-29], and also by the QCD 
sum rules in medium [[Fj. For instance, according to Brown and Rho |3l[, the masses of 


hadrons decrease at higher densities. The change of the mass leads to a change of the hadron 
densities, which should be reflected in the measured relative particle yields. The study of 
such effects is the second important issue discussed in the present paper. 


In our study of the particle multiplicities we use our own implementation of the thermal 
model, which has been constructed as a code in the MATHEMATICA language. The use of 
the symbolic manipulations allowed by the MATHEMATICA turned out to be very convenient 
in the treatment of the hadronic decays. Creation of this code was the main technical task 
connected with the present investigations. A part of the original results discussed in this 
paper has been published before in the two articles: 


I. M. Michalec, W. Florkowski, and W. Broniowski: Scaling of hadron masses and widths 
in thermal models for ultra-relativistic heavy-ion collisions , Phys. Lett. B520 (2001) 
213; |nucl-th/0103029l 


II. W. Florkowski, W. Broniowski, and M. Michalec: Thermal analysis of particle ratios 
and p± spectra at RHIC, |nucl-th/0106009| . 
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Chapter 2 

Thermal Model of Particle Production 


2.1 Historical Perspective 


The use of statistical concepts to describe particle production in hadronic collisions has a 
long history. In the early fifties Fermi [[I2j assumed that when two relativistic nucleons 
collide, the energy available in their center-of-mass system is released in a very small volume 
(due to the Lorentz contraction). Subsequently, such a dense system decays into one of many 
accessible multiparticle states. The decay probabilities were calculated by Fermi according to 
the standard rules of statistical mechanics. In the Landau’s hydrodynamic model the picture 
of the instantaneous break-up was modified by the inclusion of the expansion stage. The 
volume expansion leads to cooling and lowering of the freeze-out temperature, and causes 
that most final hadrons are the light ones. 


In the Landau hydrodynamic model [[33] the initial conditions are specified at a given 
time, when the matter is highly compressed and at rest. This description does not include 
one aspect of the particle production at high energies, namely, the fact that fast particles are 
produced later and farther away from the collision center than slow particles. This feature 
of particle production has been included in the Bjorken |34|, |35] hydrodynamic model which 
imposes boost-invariant initial conditions. 


Thermal and/or statistical interpretation of the particle production became a common 
approach for the ultra-relativistic heavy-ion collisions [9-23]. In this case large numbers 
( multiplicities) of hadrons are created, and the statistical methods seem to be appropriate. 
Of course, the produced matter may be formed in a state which is far away from the local 
equilibrium, so the use of the simple equilibrium concepts must be grounded in more involved 
studies using, e.g., the kinetic theory of hadronic matter. 
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2.2 Thermal and Chemical Freeze-Out 


In the simplest version of the thermal model one assumes that the hadronic matter created in 
nuclear collisions forms an ideal gas. A volume expansion of the gas, caused by high internal 
pressure, leads to decoupling or freeze-out of the hadrons. This process takes place when 
the mean free path of hadrons becomes compatible with the macroscopic size of the whole 
hadronic system. After the freeze-out, the particles move freely and their properties may be 
measured in the detectors. If the freeze-out process is fast, the momentum distributions of 
the particles do not change substantially, and the measured spectra have thermal shapes. 
Thus, by studying the hadron momentum distributions we may check the validity of the 
concept of thermalization. Moreover, in the thermal model the abundances of different 
hadron species are fixed by the values of a few thermodynamic parameters. Consequences 
of this fact can be verified by the measurements of the relative particle yields. 

More realistically, the hadronic system at the freeze-out can be viewed as a collection of 
subsystems. Each of such subsystems can be characterized by the individual (or local ) values 
of the thermodynamic parameters, and by the individual value of the collective velocity. The 
case described above corresponds to the situation when all local thermodynamic parameters 
are the same, and the collective velocities are small. In this approximation the measurements 
of the momentum distributions and of the relative yields should reveal the same value of 
the temperature and the same values of the chemical potentials. On the other hand, if the 
collective velocities are not negligible, the momentum spectra of hadrons are modified by 
a superposition of ’’redshift” and ’’blueshift” effects. The spectra do not reveal the true 
local freeze-out conditions anymore. It can been shown, however, that the measurement 
of the relative particle yields is not affected by the collective velocities, provided the local 
thermodynamic parameters are the same [ PU| . |5Tfl . This crucial observation makes grounds 
for a large interest in the studies of the particle ratios. In this paper we follow this trend 
and perform thermal analysis of the particle ratios measured at BNL AGS, CERN SPS, and 
BNL RHIC. 

Since the thermal fits to the particle ratios give quite large values of the optimal tem¬ 
perature, one finds typically T ~ 170 MeV, a concept of the two different freeze-outs has 
been introduced: at first the chemical freeze-out takes place, and only later the thermal (or 
kinetic ) freeze-out happens. At the chemical freeze-out the chemical contents of the hadronic 
system is established. Later, only elastic processes are possible (dominate), leading to fur¬ 
ther expansion of the system and cooling. Finally, when the system becomes sufficiently 
diluted, the true thermal freeze-out takes place (as described above). The concept of the 
two freeze-outs helps to understand the large difference between the temperature inferred 
from the investigation of the particle ratios, and the temperature inferred from the study of 
the momentum distributions. At SPS energies, the temperature of the thermal freeze-out is 
around T ~130 MeV. A detailed study of the freeze-out mechanism can be performed only 
in the microscopic framework, such as the relativistic kinetic theory. In our approach, we 
simply adopt the definition of the chemical freeze-out and check if the particle multiplicities 
measured in the ultra-relativistic heavy-ion collisions are consistent with this idea. 
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It should be noticed, however, that there exist models of hadron production in ultra- 
relativistic heavy-ion collisions, which assume that the chemical freeze-out coincides with 
the thermal freeze-out. An example of such a model is the sudden hadronization scenario 
of Ref. [^2|]. Another model with a single freeze-out is defined in Ref. |23|]. We want to 
emphasize that in this case (i.e., in the situations when a single freeze-out is considered) 
our calculations are also useful for the determination of the thermodynamic properties of 
hadronic matter. For example, the results of our analysis of the particle ratios at RHIC 
2TJ were combined with the hydrodynamic expansion in order to calculate the transverse- 


momentum spectra of hadrons [|23| . This approach yields very successful fits and shows that 
the thermal analysis of the ratios may be the first step in more complex investigations of 
other observables such as elliptic flow or HBT radii. 



chemical freeze-out 
inelastic processes cease 


thermal (kinetic) freeze-out 
elastic processes stop 


Figure 2.1: In many scenarios of the evolution of the hadronic matter in ultra-relativistic heavy- 
ion collisions, the chemical freeze-out precedes the thermal (kinetic) freeze-out. At the chemical 
freeze-out the inelastic processes stop and the chemical content of the fireball is established. Later 
on, only elastic processes are present. They lead to further expansion and cooling of the system. 
At the thermal freeze-out the elastic interactions cease, and the particles stream freely to detectors. 
A high temperature of the chemical freeze-out suggests that it happens just after hadronization of 
the quark-gluon plasma possibly created in the collisions. 


2.3 Basics of Thermal Analysis 
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2.3.1 Static Fireball 


In this Section we give the basic formulas used in the thermal analysis of the particle ratios. 
We start with the presentation of the simplest approach and assume that the (chemical) 
freeze-out takes place in a static volume. In this case the multiplicities of hadrons of species 
i are obtained from the ideal-gas expression 

N t = Vn t = Vc h J d 3 p f, (p). (2.1) 


Here V is the volume of the hadronic system at the freeze-out, Qi = 2J, + 1 is the spin 
degeneracy factor of the ith hadron, and fi(p) is the momentum distribution function. In 
the thermodynamic equilibrium the distribution functions f t (p) have a form (H = 1) 


fi ip) = 


1 


(2vr) 


exp 


' l-'ii'p) ~ Pi" 
T 


+ e 


-l 


( 2 . 2 ) 


where 

Eiip) = \/p 2 + mj (2.3) 

is the energy, T is the temperature, and p t is the chemical potential. The quantity e equals 
+ 1 for fermions (in this case ( |2.2| ) becomes the Fermi-Dirac distribution) and -1 for bosons 
(in this case (|2.2| ) becomes the Bose-Einstein distribution). The limit e —»■ 0 corresponds 
to the classical (Boltzmann) statistics. The chemical potential pi in Eq. ( |2.2| ) is a linear 
combination of the baryon, strange and isospin chemical potential 


Pi — p B Bi + p s Si + p 1 Ii. (2.4) 

Here Bi, Si, and Ii are the baryon number, strangeness, and the third component of isospin 
of the ?’th hadron, respectively. 

Introduction of the chemical potentials p B , p s and p 1 , allows us to fulfill the appropriate 
conservation laws. We assume that the strangeness of the system is zero 

E^ = 0, (2.5) 

i 

and the ratio of the electric charge to the baryon number in the hadronic fireball is the same 
as in the colliding nuclei 

EiQiNj _z 

E ? ; B,N, t A' ( } 

Here Qi is the electric charge of the ith hadron. According to the Gell-Mann - Nishijima 
formula we have 

Qi = I, + Ai±Al. (2.7) 

In the practical calculations, Eqs. (ED and m are used to fix the values of the chemical 
potentials p s and p 1 . Thus, we are left with only two independent parameters: the temper¬ 
ature and the baryon chemical potential. We note that the volume cancels in the conditions 
Q and Q, and also in the particle ratios introduced below. The values of Z and A (for 


the nuclei discussed in our paper) are given in the Table |0 
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Si 

S 

W 

Au 

Pb 

z 

14 

16 

74 

79 

82 

A 

29 

32 

184 

197 

207 

Z/A 

0.48 

0.50 

0.40 

0.40 

0.40 


Table 2.1: The values of the atomic numbers, Z, and the mass numbers, A. for the nuclei discussed 
in the present paper. 

2.3.2 Expanding Fireball 

In a more general case, when the expansion of the system at freeze-out cannot be neglected, 


one may use the Cooper-Frye formula |38|, [h| to calculate the total yield of particles of 
species i. namely 

N i = J J dE„ P 11 9i fi (V (p ■ u f - mfj . (2.8) 

Here S /t describes the element of the freeze-out hyper-surface, and tA is the local hydrody¬ 
namic four-velocity. In the local rest frame: tA = (1, 0, 0, 0) and \J(p ■ uf — rrij = p. Both 
and tA depend on the space-time position x. In general, also T and may be defined 

locally. 


Eq. 


where 


can be rewritten in terms of the number current |36 

Ni = j di:^x)j^(x), 

3i (x) = 2 f d A p9 (p°) 8 (p 2 - mfj p^ gi fi (V(p ■ uf - rri: 
d A p 6 (p°) 8 [p 2 — m 2 ') pP gi 


(2.9) 


= 2 


exp 


p ■ u — /i, 
T 


i -l 


+ e 


( 2 . 10 ) 


Eq. ( p.lOj ) is written in a manifestly Lorentz-covariant way. The step function 9(x) is defined 
by the conditions 

9{x) = 1 for x > 0, and 9{x) = 0 for x < 0. (2-11) 

In local thermal equilibrium the number current is proportional to the four-velocity, 

ft {x) = pi(x) v? (x), (2.12) 

and 


Pi (z) = (x) jf (x) 

= 2 f d A p0 (p°) 5 ( p 2 - mfj p M pi fi (V (p ■ uf - mfj 
= Pi f d 3 p fi (p) = m (T (x) , n l (x)) . (2.13) 
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a) 


u(Xj) 


• u(x 3 ) 

b) 


u(Xj) 

T 


u(x 2 ) 


' u(x 3 ) 


u(x 2 ) 


Figure 2.2: Illustration of the insensitivity of the particle ratios to hydrodynamic flow. An ex¬ 
panding system (a) can be viewed as a collection of subsystems (b) which move with collective 
velocities u(x). If the thermodynamic potentials are the same in each subsystem, the particles are 
emitted in identical proportions. In this case the measurement of the yields in the full phase-space 
reveals the local thermodynamic conditions. 

Here rii(x ) denotes the equilibrium particle density at the temperature T(x) and the chemical 
potential Hi(x). The total particle yield of species i is therefore 

Ni = j dTj^rii (T (x ), /q (x)) (x). (2.14) 

We observe that the particle ratios do not depend on a particular shape of the freeze-out 
surface as long as the local thermodynamic parameters are independent of x. In this case 
we have 

Nj_ _ Wj(T, fjJJc (s) _ m(T, ft) 

Nj nj(T, jxj) f (x) n/ij)' 

so the ratios are the same as those in a static fireball. 

The case considered above should be confronted with the analysis of the particle spectra. 
The latter are strongly influenced by the hydrodynamic flow, since the flow changes the 
individual velocities of the particles. As a consequence, the measurement of so-called inverse 
slopes of the p± spectra does not bring the direct information about the real temperature of 
the emitting hadronic source. For example, the pieces of hadronic matter moving towards 
the observer (large transverse flow at zero rapidity) seem to be hotter than those being at 
rest (small transverse flow). This is a kind of the blueshift effect, which effectively raises the 
measured inverse slope. 


16 



Strictly speaking, the arguments presented in this Section show that the particle ratios 
are insensitive to the hydrodynamic flow if all rapidities and transverse momenta of the 
particles are measured (in Eqs. ( |2.8| ) and ( |2.1G| ) one integrates over all momenta p) . From the 
experimental point of view, this requires the full 47 t acceptance of the detectors. Nevertheless, 
this condition may be relaxed for boost-invariant systems. In this case dNi/dy is independent 


of y and we have [23 


dNj/dy _ f dy dNi/dy _ Nj 
dNj/dy f dy dNj/dy Nj 


Here 


dNf 


= / 


dNi 


(2.16) 


(2.17) 


dy J d 2 p± dy 

Thus, for the boost-invariant systems it is enough to measure all particles at midrapidity, 
i.e., for y & 0. The measurements at midrapidity are also sufficient for the systems which 
are approximately boost-invariant [^3j[. This property will be used in our analysis of the 
particle ratios measured in Au + Au collisions at RHIC. 


2.3.3 Decays of Resonances 

The hadronic fireball consists of stable hadrons (with respect to strong interactions) and 
all hadronic resonances. In our analysis we use the newest edition of the review of particle 
physics [401. Practically, we include all light-flavor hadrons, i.e., hadrons containing u, d and 
s quarks. A few hadrons are not taken into account, since their properties are not known 
sufficiently well. All together we include 164 baryonic states and 110 mesonic states (treating 
separately different isospin states). When calculating the relative yields of the measured 
hadrons we should include all decay channels. This may be represented schematically by the 
expression 

p _ Nj + J2k b(k —> i) N k + YjH b (k —> /) b (/ —> i) N k + ... 

~ Nj + Ekb(k^j)N k + Ekib(k^l)b(l^j)N k + .N ’ 

where the sum over k and / includes all resonances, and b (m —> n) is the branching ratio for 
the decay process m —> n. If the decay process rn —> n does not take place, the branching 
ratio b (m —> n) is taken to be zero. The inclusion of the resonances is a very important effect, 
since the contributions from the decays are quite substantial, especially at large temperatures 
and densities. Eqs. (|2.18|) may be also used in a slightly different form, namely 


R = 


n. 


+ J2kb(k —> i)n k + 'Em b (k —> 1) b (Z —► i) n k + ... 


nj + E kb(k ^ j ) n k + E ki b{k-*l)b{l -> j ) n k + ... ’ 


where 


rii = 


p 2 dp 


(2.19) 


( 2 . 20 ) 


2vr 2 Jo exp [(E t - p*\ em Bi - p s chem Si - p{ hem Ii) /T chem ] + e ' 

To stress the fact that we calculate the ratios at the chemical freeze-out, we have marked 
the temperature and the chemical potentials with an additional label (in this way we also 
make a connection with the notation used in Refs. |42|] and |43|). 
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7T 



7t 


Figure 2.3: At the chemical freeze-out, the hadron gas consists of stable hadrons (with respect to 
strong interacions) and all hadronic resonances. The observed multiplicities obtain contributions 
from the primary particles (present in the fireball), and secondary particles coming from the strong 
decays of the resonances. In addition, the electro-weak decays should be reconstructed. In this way 
the yields of strange particles can be obtained. 


i) Isospin Symmetry 


In our approach we take into consideration all charged (isospin) states of hadrons separately. 
The non-zero value of the fitted chemical potential /i 1 leads to the splitting between the 
multiplicities of the hadrons belonging to a given isospin multiplet. Consequently, the ratios 
of such hadrons may be included in our analysis. For example, we take into consideration 
the 7r + /7r _ ratio measured at SPS and RHIC. The typical values of p 7 are of the order of 
10 MeV, thus the inclusion of the isospin chemical potential represents an important effect 
(only at RHIC p 7 is practically zero). 


A correct treatment of the different charged states requires that the branching ratios 
supplied by the Review of the Particle Physics [40] must be supplemented by the calculation 
of the probability of the transitions between such states. Since the strong interactions are 
invariant under rotations in the isospin space, the matrix elements needed for the two-body 
strong decays may be obtained directly from the Clebsch-Gordan coefficients. [] For example, 
we can consider the strong A —> 7 :N decay. Practically, this is the only decay channel of 
A, so its branching ratio is 100%. In order to include different charged states of A, 7r and 

Ror the three-body decays an additional averaging is needed, unless the full matrix element describing 
the decay is known. 


18 










the nucleon, N = ( p,n ), we combine the isospin states 1 and 1/2 into the isospin states 3/2: 


A ++ 

A + 

A 0 

A~ 


1 1 1 .3 3 \ , , , , 

’2 2 2 2 / 71 +P) = 7T +P ’ 

{ n+ + n ) + ( 10 ’ll\ll) (*°+p) 

i (tt + + nj + jj (tt° + p) , 

(l . 1 ,- | |_i )(# . +J() + (l0t i_i | |_i)(^ B ) 

^ (tt _ + p) + | (*° + n) , 

A , 1 1 ,3 3\ 

\ ’2 2 '2 2 / y J 


ii) Experimental Branching Ratios 


The case A —> nN is very comfortable, since only a single decay channel is present and 
the branching ratio is obviously well known. In majority of the cases we deal with, different 
decay channels appear and their properties (branching ratios) are sometimes not well known. 
As a rule we disregard all decays with the branching ratios smaller than 1%. In addition, if 
the decay channels are described as dominant, large, seen, or possibly seen, we always take 
into account the most important channel. If two or more channels are described as equally 
important, we take all of them with the same weight. For example /o(980) decays into tttt 
(according to ||(J this is the dominant channel) and KK (according to [|dCj|] this is the seen 


channel). In our approach, according to the rules stated above we include only the process 
/o(980) —> 7T7T. Similarly, ao(1450) has three decay channels: rjTT (seen), 7tt/(958) (seen), 
and KK (again seen). In this case we include all three decay channels with the weight 
(branching ratio) 1/3. 

Another difficulty is that the branching ratios are not given exactly (instead of one value, 
the whole range of acceptable values is given) and the sum of the branching ratios may differ 
significantly from 1. In this case we take the mean values of the branching ratios. Since 
we require that their sum is properly normalized, sometimes we are forced to rescale all the 
mean values in such a way that their sum is indeed 1. To illustrate this problem we present 


our analysis of the decays of the A(1600) resonance in Table |2[^. Below, as an example we 
give the final branching ratios of the A(1600) ++ resonance: 


A(1600) 


++ 


0.175 (vr + +p) 

+ 0.504 [- (A + 

. 5 v 

+ 0.115 (p + p + ) 


+ ) + 


(a ++ 


+ 7r 


0.206 


A(1400) + + 7 
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A(1600) decay modes 

branching ratio 

averaged 
branching ratio 

rescaled 
branching ratio 

Nn 

10-25 % 

17.5 % 

17.5 % 

Nnn 

75-90 % 

82.5 % 

82.5 % 

An 

40-70 % 

55.0 % 

50.4 % 

Np 

< 25 % 

12.5 % 

11.5 % 

A(1440)tt 

10-35 % 

22.5 % 

20.6 % 

Ay 

0.001-0.02 % 

— 

— 

Ay, hclicity=l/2 

0.0-0.02 % 

— 

— 

Ay, hclicity=3/2 

0.001-0.005 % 

— 

— 


Table 2.2: The branching ratios for the decays of A(1600). The experimental information (second 
column) is averaged (third column) and rescaled (fourth column) in order to achieve the correct 
normalization. Similar procedure is also applied for other decays whose branching ratios are poorly 
known. 

We note that the correct normalization of the branching ratios is crucial for the fulfillment of 
the conservation laws, and consequently for the determination of the optimal thermodynamic 
parameters. 



Figure 2.4: A schematic view of the weak decays of strange baryons. 


iii) Weak Decays 

An important part of the thermal analysis is also the correct treatment of the contributions 
from the electro-weak decays. For example, the final (measured) multiplicity of pions is 
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modified by the decays of A, E and 5. Since the abundances of A, E and S in a hot and 
dense matter are not negligible, their decays modify the final pion multiplicity. To study 
such effects, we do our calculations usually in three different ways. In the first version we 
include all contributions from the weak decays. In the second version we assume that there 
is no feedback from the weak decays - this case corresponds to the ideal situation when all 
weak-decay channels can be experimentally disentangled. Finally, in the third version we 
assume that the feedback from the weak decays is at the level of 50%. Here we follow the 
concept of Ref. 


13 


To be more precise, in our description of the weak decays the labels a, b, c and d are 
introduced (see Fig. |2.4| ). The label a specifies the particle multiplicity which includes the 
two contributions: the first one from the primordial particles, and the second one from all 
strong-decay channels. At the a-level our system includes: 7r, rj, K, N, A, E, 5 and fh The 
label c characterizes the abundances corrected for the and E decays, whereas label d means 
that the decays of 5 are also included. At the d-level we deal with: n, r /, K, N and A. The 
final abundances of pions, kaons and nucleons (including all electro-weak decays) are denoted 
by the label b. 


iv) Resonance Mass Spectrum 


In Refs. |45| , |46| , |47j the arguments have been presented that the mass spectra of baryons 
and mesons behave differently: the baryon spectrum grows more rapidly than the meson 
spectrum. As a consequence, the Hagedorn temperature (i.e., a scale describing the expo¬ 
nential growth of the spectrum) is different for baryons and mesons. With the use of a 
simple exponential formula, the spectra of baryons and mesons may be well parameterized 
as follows: 

777 / 

PB,Ai{m) = A BM ex p 


rpB,M 

L 1 H 


( 2 . 21 ) 


with the baryon Hagedorn temperature 


T| = 186 MeV, 


( 2 . 22 ) 


the meson Hagedorn temperature 


— 311 MeV, 


(2.23) 


and the normalization constants: Ab = 4.41 GeV and Am = 0.11 GeV [B3. 


2.3.4 Sketch of x 2 Method 

The optimal values of the temperature T chem and the baryon chemical potential /i^ em are 
fitted by minimizing the expression 
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(2.24) 


2 ^ (. RT P - R{ herm ) 2 

x h °i 

where R e ^ p is the kth measured ratio, ay is the corresponding error, and R t ^' er,n is the same 
ratio determined from the thermal model. The total number of different ratios included in 
the analysis is denoted by n. 

Introducing the short-hand notation ay = T c h em and ct 2 = Pchem we ma y write 


X (®) X (®min) T 


CTn 


) V (ft C^rnin) ) 


(2.25) 


where a = (ap,^), «mm is the optimal pair of the parameters, and V is the variance 
matrix of the parameters a. We note that y 2 (ct) has a minimum at a = a m i n , so the first 
derivatives of y 2 (a) vanish at this point. If F (a) is some function of the fitted parameters 
a, the variance of F is given by 


1 j hdOLmdou mn • 


(2.26) 


In the special cases F — ay and F = a 2 , Eq. (|2.26D can be used to make an estimate of the 
errors of the fitted temperature and baryon chemical potential 


(ATchemY — Ell, 




l chem J 


= Vo,. 


(2.27) 


Equation (|2.27|) will be used in our analysis with the coefficients Vii and V 22 determined 
numerically. 


2.3.5 Thermodynamic Characteristics of Freeze-Out 

It is very interesting to know the thermodynamic properties of the hadronic matter at freeze- 
out. Of the particular interest is the energy density of the hadronic system. Its closeness to 
the critical energy density for the deconfinement phase transition (found by the Monte-Carlo 
simulation of QCD on the lattice |5D|]) may indicate that such a phase transition indeed 
took place in the ultra-relativistic heavy-ion collisions. In our approach the knowledge of the 
temperature and the chemical potentials allows us to calculate all intensive thermodynamic 
quantities. In particular, we calculate the energy density 


= J2 £ i = J29i / d 3 P Ei (p) fi (p) , 


(2.28) 


pressure 


p = J2 p i = J2^ d 3 P 


3 ^ p 


3 Ei (p) 


fi (p), 


(2.29) 
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and the baryon, strangeness, and isospin densities 


p B = J2di J d 3 pBifiip ), 

i 


(2.30) 

Ps = J2di J d 3 pSifi(p ), 

i 


(2.31) 

Pi = J2di J d 3 plifi(p). 


(2.32) 

The entropy density is determined from the Gibbs identity 



s + P — p B p B ~ PsPs ~ Pi pi 

S ~ T 


(2.33) 

It is interesting to discuss the case of the classical statistics separately (this is achieved 
by taking the limit e —> 0 in the equilibrium distribution functions (^.2|)). In this case the 
thermodynamic quantities can be expressed in terms of the modified Bessel functions 

7- / 1 2 nn! —n f°° j (2 2\ n ~ 1 /2 _ T 

Kix)= (2n)\ X id' -) e • 


(2.34) 

In particular, the particle densities can be found from expression 



n i - 27r2 e M ' /T mfrW 2 ( T ) , 


(2.35) 

and the corresponding energy density and the pressure are determined by the following two 
equations 

+ - 2 V‘ /Tm?TA ' 3 (r) 


(2.36) 

P ‘ ~ 


(2.37) 


One can notice that Eqs. (|2.35|) and (|2.37|) yield the equation of state of the relativistic ideal 
gas of classical particles (k B = 1) 

P = nT. (2.38) 


We note that it is the same as the equation of state of the non-relativistic ideal gas. 


It has been pointed out by Cleymans and Redlich [lj], that for a variety of different 
colliding systems the average energy per hadron at the chemical freeze-out is very close to 1 
GeV 

£ Ei Si 




1 GeV. 


(2.39) 


n Ei rii 

In our calculations we check this relation. In addition, we calculate separately the average 
energy of baryons and mesons defined as 


r B = 


SB 

n B 


E Si 

baryons 

E n i 

baryons 


r M = 


Sm 

Um 


E Si 

mesons 

E 

mesons 


(2.40) 
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2.4 Modifications of Thermal Analysis 

2.4.1 Finite-Size Corrections 


The particle densities defined by the ideal-gas expression correspond to the thermodynamic 
limit, i.e., they are calculated in the limit V —> oo, N —> oo, N/V = const. In realistic 


situations, the hadronic systems have finite volumes and the approach based on Eq. (|2.20|) 
may be not sufficiently accurate. To include the effects of the restricted volume, one can 
use the modified versions of Eq. ( |2.20| ) . In practice, such modifications are known only for 
simple geometries. In the special case of the spherical symmetry, when the hadronic system 
forms a fireball of radius R s , a correction of Eq. ( |2.20|) may be achieved by the change of 
the level density |5D 

3vr 1 


dp —> dp 1 — 


4 pR s + (pRsf 


+ 


Hence, the particle densities may be calculated from the formula 


(2.41) 


m = 


9i 

2vr 2 


dp p 2 


3np 1 

4iF + m + "\ 


exp 


Pi Pi,chem^ _ ^ 


- chem 


(2.42) 


Clearly, for very large radii, R s —> oo, Eq. (|2.42| ) is reduced to the standard formula ( |2.20|) . 
The practical calculations B0 indicate that the finite-size corrections affect mainly the 
absolute densities. The particle ratios are not changed much due to this effect. 


2.4.2 Excluded-Volume (van der Waals) Corrections 


The excluded-volume corrections account for the finite volumes of hadrons. In a very dense 
hadronic matter the particles with finite size may overlap. This leads to strong repulsive 
forces and the simple approach based on Eqs. ( |2Tp and ( |2.2|) should be modified (by 
definition Eqs. (|2.1| ) and (|2.2|) describe non-interacting particles). Of course, the overlapping 
of hadrons may be the first step towards the deconfinement phase transition and creation of 
the quark-gluon plasma. In any case, however, the repulsive part of the nuclear force at short 
distances should be included in the realistic description of the hadron thermodynamics. 


A fully consistent (from the thermodynamic point of view) method to include the ex¬ 
cluded volume corrections was introduced by Yen, Gorenstein, Greiner and Yang [^T|. In 
their approach one calculates the modified pressure P defined by 


P(T,pi,p2, •••) — P(T, pi,p2, ■■■), (2.43) 

where P is the pressure of the ideal gas, as defined by Eq. (|2.29|) , and /2i are the modified 
chemical potentials 
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Pi = Pi~ViP{T,p i,// 2 , ■■■)■ (2.44) 

Iii Eq. (|2.44|) i>i is the particle eigenvolume |5T|. We note that Eq. (|2.43|) is a non-linear 


equation for P, which can be solved by the iterative method. The particle densities fq are 
the derivatives of P with respect to the chemical potentials p t 


dP ^ dP d/j,j _ 

n ^ = — =l^ — = zL n i ~ v i n i> ’ 




djlj dfii 


(2.45) 


3 


which gives 


rii(T, Hi) 


rti(T , /Xj) 

1 + J2 j v j n j (T, Jij)' 


(2.46) 


We observe that the modified densities are smaller than the initial densities. The reasons 
for this are twofold: firstly, the denominator appearing on the right-hand side of Eq. ( |2.46|) 
is always larger than one; secondly, the densities n t are calculated at smaller values of the 
chemical potential, since we have always pi < /q. 


For our investigation of the relative particle yields, it is important to realize that the 
denominator in Eq. ( |2.46|) is the same for all hadron species. Hence, it cancels in the 
particle ratios. In addition, for the classical (Boltzmann) statistics and identical particle 
volumes of all hadrons (n* = v), the modifications of the chemical potential factorize, and 
again cancel in the ratios. Thus, in most cases the excluded volume corrections do not affect 
the thermal analysis of the particle ratios and have no impact on the fitted value of the 
optimal temperature and baryon chemical potential 


'W'i(T, pf) ^ ?lj(P,/lj) 

nj(T, pf) ni(T, pf) 

Note, however, that the actual particle densities should be calculated from Eq. 
the modified chemical potential Jii. 


(2.47) 
(|2.44j ) with 


The case of the Boltzmann statistics is also interesting, since it leads directly to the van 
der Waals equation of state. To see this feature, we first write 


P = exp 



(2.48) 


A derivative of Eq. ( |2.48| ) with respect to the chemical potentials yields 


P (1 — vn) = nT , n — y~) rq, 

i 


(2.49) 


which is a typical excluded volume correction. 
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2.4.3 Corrections for Undersaturation of Strangeness 


The chemical equilibrium between strange and non-strange particles is more difficult to 
achieve because the production of ss pairs proceeds usually at a slower rate than the pro¬ 
duction of uu and dd pairs. Moreover, there are no strange quarks at the beginning of 
the collision process. In Ref. |12| Rafelski has introduced an extra parameter to account 


for possible undersaturation of strangeness. In that approach one calculates the hadronic 
densities from the formula 


tii = 


9i 

27r 2 


dpp 2 


7~«exp( Ei ^ ch ™ ) +e 


T, 


chem 


■ 


-i 


(2.50) 


where 7 is the strangeness saturation factor (0 < 7 < 1) and s t is the number of valence 
strange quarks in the ith hadron. 


2.5 Thermal Model in this Work 


In the following Sections we apply the thermal model to describe the particle ratios measured 
in different experiments with ultra-relativistic heavy-ions. We use the simplest version of 
the model, and usually neglect the corrections discussed in Sections [2.4.1| - |2.4..‘j| . In other 
words, we assume the full chemical and thermal equilibrium, and use the ideal-gas formulas. 
The role of the corrections defined in Sections [2.4. 1| - [2.4.3| is discussed in special cases. Our 
method of dealing with the resonance decays was presented in detail in Section [2.3.3| . The 
results of our calculations are denoted in the text by the label TM. 
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Chapter 3 

Reactions with Si and S Beams 


In this Chapter we use our implementation of the thermal model to analyze some of the first 
experiments with the ultra-relativistic heavy-ions: the experiments with the 28 Si beams at 
BNL AGS, and the experiments with 32 S beams at CERN SPS. The Alternating Gradient 
Synchrotron at BNL accelerated 28 Si ions to the energy of 14.6 GeV per nucleon (in 1986), 
whereas the Super Proton Synchrotron at CERN accelerated 32 S to the energy of 200 GeV 
per nucleon (in 1987 and later in 1990). For many years the results from sulphur-induced 
reactions were the main source of the data on hadron production at true ultra-relativistic 
energies. 


3.1 Si + Au(Pb) Collisions at BNL AGS 


The heavy-ion collisions with Si beams were analyzed in a thermal model by Braun-Munzinger, 
Stachcl, Wessels, and Xu 0 (for Au and Pb targets), and later by Cleymans, Elliot, Satz, 
and Thews [|TD[ (for Au targets). In this Section we are going to compare their results with 
our analysis of the particle ratios. 


In Ref. H the finite-size and the excluded-volume corrections were included, and the par¬ 
ticle ratios were calculated for two different temperatures, T chem = 120 MeV and T chem = 140 
MeV, and the fixed baryon chemical potential ^hem = 540 MeV. The range of temperatures 
considered in [j9j was motivated by the experimental spectrum of the A(1232) resonance. 
The value of the baryon chemical potential was constrained by the measured pion to nucleon 
ratio. For a given temperature and baryon chemical potential the strangeness chemical po¬ 
tential was adjusted to give zero net strangeness (compare Eq. ([231) ). On the other hand, 
the isospin chemical potential was neglected. The results of Ref. |9j are shown in Table E3 
in the fifth (BM1) and seventh (BM2) column: one can see that the overall agreement with 
the data (second column) has been achieved, but the statistical significance of this result is 
rather poor, which is indicated by the very large values of y 2 . In the two discussed cases 
one Ends: y 2 /n = 8.5 and y 2 /n = 12.2. 
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Si+Au(Pb) 

experiment 

TM 

TM 1 

BM 1 

TM 2 

BM 2 

Tchem [MeV] 


136±2 

120 

120 

140 

140 

Pchem [MeV] 


593±9 

540 

540 

540 

540 

Pchem [MeV] 


152 

109 

108 

135 

135 

Pthem [MeV] 


-14 

-11 

0 

-13 

0 

x 2 /n 


6.3 

9.3 

8.5 

22.1 

12.2 

■k/ (p + n) 

1.05(5) 

1.20 

1.36 

1.29 

1.45 

1.34 

P/p 

4.5(5)T0~ 4 

3.6T0 -4 

2.1-10 -4 

1.47-10 -4 

9.0-10 -4 

5.8T0 -4 

K + /n + 

0.19(2) 

0.23 

0.23 

0.23 

0.24 

0.27 

K~ / 7T~ 

3.5(5)T0~ 2 

3.3T0 -2 

3.9-10 -2 

5.0-10” 2 

4.3-10 -2 

6.2-10” 2 

K“/n+ 

9.7(15)T0 -2 

14-10 -2 

14-10 -2 

14-10~ 2 

15-10” 2 

16T0 -2 

K+/K~ 

4.4(4) 

5.5 

4.6 

4.6 

4.4 

4.3 

A/(;p + n) 

8.0(16)T0 -2 

7.2-10 -2 

7.6-10 -2 

9.5-10 -2 

8.6-10 -2 

11-10~ 2 

A/A 

2.0(8)T0 -3 

2.3T0 -3 

1.0-10” 3 

0.88T0 -3 

4.4-10 -3 

3.7T0 -3 

0/ (K+ + K~) 

1.34(36)-10~ 2 

2.89T0 -2 

2.4-10 -2 

2.4-10 -2 

3.5-10 -2 

3.6T0 -2 

S-/A 

0.12(2) 

0.05 

0.056 

0.064 

0.059 

0.072 


Table 3.1: Thermal analysis of the particle ratios measured in Si + Au(Pb) reactions. The 
experimental ratios are used in the form prepared by Braun-Munzinger et al. in Ref. ||. The 
optimal values of the thermodynamic parameters found in our approach are presented in the third 
column (TM). The results of Ref. p] are presented in the fifth and seventh column (denoted by 
BM1 and BM2, respectively). For comparison the results of our model for the same input values 
°f T c hem and Uchem are shown in the fourth and sixth column (TM1 and TM2, respectively). In 
our calculations (TM,TM1,TM2) we include the weak decays at the level of 50%. 


We have reanalyzed the experimental data (as compiled in ||) in our model assuming 
the full chemical equilibrium and neglecting the finite-size and excluded-volume corrections. 
The results of our fit are shown in the third column (TM). We have found a slightly better 
fit with x 2 / n = 6.3 for T chem = 136 MeV and = 593 MeV. Still, the value of x 2 / n is 
not satisfactory. As expected, using the values of Ref. || for T chem and Pchem as an input in 
our model (without the fitting procedure), we find the ratios corresponding to higher values 
of x 2 / n ■ These results are shown in the fourth and sixth column denoted by TM1 and TM2, 
respectively. We note that our determination of the strange chemical potential agrees with 
the result of || in these two cases. 


In Table [3]2| we show the values of the thermodynamic parameters which follow from our 
analysis of the particle ratios. Our best fit gives the energy density e = 0.6 GeV/fm 3 , and the 
baryon number density ps = 0.34 frn -3 (we use here Eqs. Q2.28|) and (|2.30|) ). In all considered 
cases the pressure satisfies the condition P = (hm + %) T, which indicates that the matter 
at the chemical freeze-out behaves like an ideal classical gas (compare Eq. (|2.38|) ). In this 
situation the excluded-volume corrections cancel in the particle ratios and do not affect the 



































Si+Au(Pb) 

TM 

TM 1 

TM 2 

T chem [MeV] 

136 

120 

140 

Vchem [MeV] 

593 

540 

540 

Uchem [MeV] 

152 

109 

135 

Vchem [MeV] 

-14 

-11 

-13 

Eb [GeV/fm 3 ] 

0.48 

0.14 

0.42 

Em [GeV/fm 3 ] 

0.12 

0.05 

0.14 

e [GeV/fm 3 ] 

0.60 

0.19 

0.56 

P [GeV/fm 3 

0.07 

0.02 

0.07 

p B [1/frn 3 ] 

0.34 

0.11 

0.29 

s [1/fm 3 ] 

3.4 

1.3 

3.3 

ns [1/fm 3 

0.34 

0.11 

0.29 

n M [1/fm 3 ] 

0.17 

0.09 

0.19 

r B [GeV] 

1.4 

1.3 

1.4 

r M [GeV] 

0.7 

0.6 

0.7 

r [GeV] 

1.2 

1.0 

1.1 


Table 3.2: Thermodynamic parameters at the chemical freeze-out, as inferred from the analysis 
of the particle ratios in the Si+Au(Pb) collisions at AGS. The results correspond to the Table] 
where the fitted ratios are shown. 


optimal values of the temperature and the baryon chemical potential (as discussed in the end 
of Section |2.4.2|) . The classical form of the equation of state is a consequence of the fact that 
most of the particles included in the thermal approach are heavy resonances, which are well 
described by the classical distribution functions. The important role of heavy resonances in 
the thermal analysis is reflected also by the fact that the final density of pions (resulting 
from the decays of resonances) is much higher than the primordial density. Initially, the 
pion density is 0.09 fm -3 . With the inclusion of the strong decays of the resonances, the 
pion density increases up to 0.35 fm -3 . We also note that the contribution from the decays 
of the p-mesons to the final pion density is at the level of 15%. This means that most of the 
extra pions is produced by the decays of heavier resonances. The last raw in Table shows 
the Cleymans-Redlich ratio calculated for three different choices of T c h em and /i^ em . In all 
cases we observe that r ps 1 GeV. It is interesting to notice, however, that r B is significantly 
larger than Em- This feature reflects a different behavior of the mass spectra of baryons and 
mesons, as displayed by Eqs. ( |2.21|) - ( j2.23|) . 

The large values of X 2 /n, shown in Table 371, indicate that thermal description of the 
particle yields does not work well for reactions with the silicon beams. Even the ratio of pions 
to nucleons does not come out accurately, giving quite large contribution to y 2 /n. In the 
thermal analysis of Ref. || the isospin chemical potential is zero, so the ratios such as n + /p 
or 7r /p are necessarily the same. The experimental estimate of Ref. [9] for any of these 
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Si+Au 

experiment 

TM 

CEST 

T chem [MeV] 


108±5 

110±5 

Vchem [MeV] 


540±10 

540±20 

Vchem [MeV] 


93 


^chem [MeV] 


-9 


X /n (n = 5) 


0.5 

0.6 

7T + /p 

0.80T0.08 

0.84 

0.87 

K + / 7T+ 

0.19T0.02 

0.20 

0.21 

K+/K~ 

4.40T0.40 

4.50 

4.51 

A/p 

0.20T0.04 

0.14 

0.16 

K~/ 7T- 

0.035T0.005 

0.035 

0.038 

S-/A 

(1.2±0.2)-10“ 1 

5.2-10 -2 

4.9T0 -2 

0/7r + 

(4.5±1.2)-10“ 3 

4.2-10 -3 

4.6-10" 3 

p/p 

(4.5±0.4)-10~ 4 

7.2-10 -5 

7.2-10 -5 

A/A 

(2.0±0.8)T0“ a 

3.2-10 -4 

3.4-10" 4 


Table 3.3: Thermal-model analysis of the particle ratios for Si + Au collisions. The experimental 
data are taken from the paper by Cleymans et al., see Ref. [jlO:]. The results of Ref. [10 are given 
in the last column (CEST) and compared to our calculation (TM). In the calculation of \ 2 only 
five ratios of the most abundant hadron species are included: vr + /p, K + /-rr + , K + /K ~, A/p and 
K~/n~. The remaining four ratios (the bottom four raws) are treated here as an output of the 
thermal model. In our calculation the weak decays were included again at the 50% level. 


ratios is (2/3) tt/(p + n) = (2/3) 1.05 ~ 0.7. This result was updated in Ref. |l(J, where the 


collisions with Au targets are included only, and the following values are used: tt + / p = 0.8 
and 7 t~ /p = 1.0. The approach of Ref. m has the non-zero isospin chemical potential, so 
differences in the abundances of the states belonging to the same isospin multiplets can be 
easily included. In addition, the authors of Ref. m decided to apply the thermal model to 


the particle species which are most abundant. Consequently, their fit includes only five ratios: 
7r + /p, K + / 7r + , K + /K~, A/p and K~ /7r — . Note, that the ratio tt~/ p is not independent and 
it is not included in the calculation. 


The results of the thermal fit of Ref. |T0| are presented in Table |3.3| in the last column 


denoted by CEST. The experimental data in the second column are taken also from Ref. 
eg. We have run our code and obtained the results shown in the third column denoted by 
TM. First of all, we see that the five ratios included in the analysis are quite well reproduced 
in the thermal model (x 2 /n ~ 0.5). Our calculation fully confirms the result of Ref. m 
which gives: T c h em = 110 MeV and Pchem = 540 MeV. The four remaining ratios of less 
abundant hadron species are treated in this case as an output of the thermal model. They 
are shown in the four bottom raws of Table |3.3| . Naturally, the agreement of the model in 
this case is not as good as for the five “input” ratios. 
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Si+Au 

experiment 

TM 

CEST 

T chem [MeV] 


127±3 

110±5 

Pchem [MeV] 


533±13 

540±20 

^chem [MeV] 


116 


h, [hem [MeV] 


-11 


X / n (n = 9) 


3.1 

12.2 

7 r + /p 

0.80T0.08 

0.87 

0.87 

K + /n + 

0.19T0.02 

0.23 

0.21 

K+/K~ 

4.40T0.40 

4.42 

4.51 

A Ip 

0.20T0.04 

0.17 

0.16 

K~/tt~ 

0.035T0.005 

0.042 

0.038 

S-/A 

(1.2±0.2)-10 1 

5.8-10” 2 

4.9T0 -2 

(j)/ll + 

(4.5±1.2)-10- a 

8.2-10 -3 

4.6T0 -3 

P/P 

(4.5±0.4)-10“ 4 

4.1-10" 4 

7.2-10" 5 

A/A 

(2.0±0.8)-10“ a 

1.9T0 -3 

3.4-10” 4 


Table 3.4: Thermal-model analysis of the particle ratios for Si + Au collisions. The experimental 
data are taken from jlCj . In this case, our calculation (TM) includes all 9 ratios in the construction 
of x 2 - The last column repeats the results of Cleymans et al., see Table which were obtained 
by fitting only the first 5 ratios. The value of x 2 / n is calculated in both cases (TM and CEST) for 
all 9 ratios. 


In Table T4 we show the result of the thermal analysis obtained in the case when all 
nine ratios are included in the calculation of y 2 . Our calculation gives T c h em = 127 MeV and 
l^chem = 533 MeV. As expected, the quality of the fit, x 2 / n ~ 3, is worse than that obtained 
in the previous case with only five ratios included. On the other hand, it is better than the 
quality of the fits presented in Table |3.1| for Si + Au(Pb) collisions. The better agreement 
is caused by the use of the updated ratios of pions to nucleons in Ref. 
consistent with the thermal picture. 


101, which are more 


3.2 S -(- Au(W,Pb) Collisions at CERN SPS 


The thermal analysis of the S + Au(W,Pb) collisions was performed by Braun-Munzinger, 
Stachel, Wessels, and Xu in Ref. [[□]]. In Table [15] we show their results (BM1, BM2) 
together with our analysis based on the assumption of the full chemical and thermal equi¬ 
librium (TM, TM1, TM2). In our approach we have again neglected the finite-size and the 
excluded-volume corrections. The presentation of the results in Table [h5] is analogous to 
that given in Table [1.1|. We see again that the values of y 2 /n are quite large, indicating that 


the thermal model describes the measured ratios only in a qualitative way. 
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S + Au(W,Pb) 

experiment 

TM 

TM 1 

BM 1 

TM 2 

BM 2 

Tchem [MeV] 


179±2 

160 

160 

170 

170 

Pchem [MeV] 


199±4 

170 

170 

180 

180 

Pchem [MeV] 


55 

37 

38 

45 

47 

Rcfcem [MeV] 


-7 

-4 

0 

-5 

0 

X / n 


4.2 

11.3 

7.7 

5.6 

6.6 

K"JA d 

0.88(10) 

1.29 

1.68 

1.57 

1.46 

1.36 

K a ,J A„ 

4.6(10) 

5.4 

7.8 

7.3 

6.15 

5.7 

Ac/A c 

0.20(1) 

0.21 

0.19 

0.20 

0.21 

0.23 

A d/ A d 

0.19(4) 

0.24 

0.21 

0.22 

0.24 

0.24 

“c/Ac 

0.095(6) 

0.112 

0.12 

0.12 

0.16 

0.12 

S+/A c 

0.21(2) 

0.20 

0.18 

0.20 

0.19 

0.21 

p/lt + 

0.18(3) 

0.29 

0.22 

0.17 

0.25 

0.19 

(. h + - h~)/h~ 

0.15(1) 

0.19 

0.14 

0.18 

0.16 

0.21 

(p - P) / h ~ 

0.15(2) 

0.20 

0.15 

0.13 

0.17 

0.14 

(. h + - h~) / (h+ + h~) 

0.088(7) 

0.085 

0.067 

0.084 

0.075 

0.094 

P/p 

0.12(2) 

0.14 

0.14 

0.13 

0.15 

0.14 

p/n- 

0.024(9) 

0.037 

0.029 

0.022 

0.035 

0.027 

T}/ 7T° 

0.15(2) 

0.11 

0.12 

0.12 

0.11 

0.12 

H ip + w) 

0.080(20) 

0.055 

0.053 

0.11 

0.055 

0.12 

(K* + K~)/(2K») 

1.07(3) 

1.03 

1.03 

1.05 

1.03 

1.06 

K+/K~ 

1.67(15) 

1.58 

1.44 

1.46 

1.49 

1.53 

K! A 

1.4(1) 

1.6 

2.12 

1.74 

1.84 

1.50 

Kg / A 

6.4(4) 

7.1 

10.3 

8.5 

8.1 

6.6 

A /(p - p ) 

0.45(4) 

0.43 

0.44 

0.67 

0.45 

0.69 

A/p 

0.80(30) 

0.62 

0.55 

0.38 

0.59 

0.41 

A/A 

0.207(12) 

0.226 

0.206 

0.20 

0.226 

0.23 

S“/A 

0.066(13) 

0.114 

0.120 

0.12 

0.118 

0.12 

S+/A 

0.127(22) 

0.185 

0.172 

0.20 

0.178 

0.21 

E + /E~ (WA85) 

0.45(5) 

0.37 

0.30 

0.31 

0.34 

0.36 

E+/E~ (NA36) 

0.276(108) 

0.37 

0.30 

0.31 

0.34 

0.36 

(0++Q-) 

(^ + +^-)_ 

0.8(4) 

0.2 

0.2 

0.17 

0.18 

0.19 


Table 3.5: Thermal analysis of the particle ratios measured in S + Au(W,Pb) reactions. The 
experimental ratios are taken from Ref. 0 - The optimal values of the thermodynamic parameters 
found in our approach are presented in the third column (TM). The results of Ref. [11] are presented 
in the fifth and seventh column (denoted by BM1 and BM2, respectively). For comparison the 
results of our model for the same input values of T c h em and Pchem are shown in the fourth and 
sixth column (TM1 and TM2, respectively). Except for the first six measured ratios (where the 
weak decays have been disentangled), we include the weak decays at the level of 50%. 
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S + Au(W,Pb) 

TM 

TM 1 

TM 2 

T chem [MeV] 

179±2 

160 

170 

Vchem [MeV] 

199±4 

170 

180 

hchem [ MeV ] 

55 

37 

45 

Vchem [MeV] 

-7 

-4 

-5 

e b [GeV/fm 3 ] 

0.44 

0.15 

0.26 

Em [GeV/fm 3 ] 

0.64 

0.30 

0.45 

e [GeV/fm 3 ] 

1.1 

0.45 

0.71 

P [GeV/fm 3 ] 

0.16 

0.07 

0.11 

p B [1/frn 3 ] 

0.20 

0.07 

0.12 

s [1/fm 3 ] 

6.7 

3.1 

4.7 

ns [1/fm 3 

0.27 

0.09 

0.16 

n M [1/fm 3 ] 

0.65 

0.35 

0.49 

r-B [GeV] 

1.6 

1.6 

1.6 

r M [GeV] 

1.0 

0.8 

0.9 

r [GeV] 

1.2 

1.0 

1.1 


Table 3.6: Thermodynamic parameters at the chemical freeze-out, as inferred from the analysis 
of the particle ratios in the S+Au(W,Pb) collisions at SPS. The results correspond to the cases 
described in Table jP] . 

In Table |3.6| we give the values of the thermodynamic parameters at the freeze-out. 
They exhibit similar features to those found in the case of the Si + Au collisions discussed 
in the previous Section. In particular, we observe the classical equation of state, quite large 
energy density and baryon number density, different Cleymans-Rcdlich ratios for mesons and 
baryons. 


3.3 Concluding Remarks 

Thermal analysis of the particle ratios measured in the collisions of relatively lighter nuclei, 
such as Si or S, leads only to a qualitative agreement with the data. This may indicate 
that the produced systems are not in the full chemical equilibrium, or they are too small 
for the thermodynamic concepts to be applicable. It is also possible that the errors of the 
experimentally determined ratios are underestimated (the true systematic errors may be 
larger, note the discrepancy seen in the measurement of the S + /S~ ratio done by different 
groups). Our calculations based on the y 2 method agree well with the results of Ref. 
where only a limited number of the particle ratios was studied. On the other hand, we 
observe appreciable differences between our fits and the global parametrization of the data 
proposed in Refs. ! El¬ 


io 
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Chapter 4 

Pb + Pb Collisions at CERN SPS 


Lead on lead collisions at 158 GeV per nucleon were the first really heavy -ion collisions at 
fully relativistic energies. In this case we dealt, for the first time, with large volumes and 
life-times of the reaction region. Simple estimates based on the Bjorken hydrodynamic model 
[ |34| indicate that matter with the energy density of a few GeV/fm 3 was created at the early 
stage of the central collisions. Such energy density exceeds the critical energy density for the 
phase transition from a hadron gas to a quark-gluon plasma |49] 
of the plasma could occur in these reactions |]52] . 


so the transient existence 


In the first part of this Chapter we show the results of the thermal analysis of the particle 
ratios, which is based on our own implementation of the thermal model. Subsequently, we 
discuss the impact of the possible in-medium modifications of hadron masses and widths on 
the thermal fits. 


4.1 Thermal Analysis of Particle Ratios 


With the same implementation of the thermal model as that used in the previous Chapter 
to describe the collisions of lighter nuclei, we have studied the particle ratios measured for 
Pb + Pb collisions at 158 GeV per nucleon. We have used the data in the form compiled 
by Braun-Munzinger, Heppe, and Stachel in Ref. |[P§ . The authors of Ref. |L3| argued 
that the particle ratios can be well described in the framework of the thermal model. Our 
calculations support this point of view - our fits to the particle ratios are shown in Table | 


To check the validity of our results, we did our calculations in three different ways, treating 
differently contributions from the weak decays. The three methods were introduced in the 
end of Section (see also Fig. [2.4| where a scheme of the weak decays is presented). 


The first 6 ratios listed in Table ^4] were measured in such a way that the feedback from 
the weak decays was well established. Consequently, these ratios enter in the same way the 
three versions of our calculation. The remaining 14 ratios are treated differently in each 
case. For example, the A jh~ ratio equals: A d/h^ in the first version (maximal feeding, 
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Pb+Pb 

experiment 

max. feeding 

no feeding 

50% feeding 

Tchem [MeV] 


164+3 

170+3 

168+3 

/C„, [MeV] 


234+7 

250+8 

244+8 

Khem t MeV ] 


56 

65 

62 

Kchem [ Me V] 


-8 

-9 

-9 

X /n 


1.4 

1.8 

0.9 

(.Pa - Pa) IK 

0.228(29) 

0.209 

0.233 

0.224 

PalPa 

0.055(10) 

0.062 

0.058 

0.060 

Pb/Pb 

0.085(8) 

0.083 

0.080 

0.081 

Ac/Ac 

0.131(17) 

0.118 

0.118 

0.118 

“c/Ac 

0.110(10) 

0.107 

0.104 

0.105 

S+/A c 

0.206(40) 

0.202 

0.210 

0.207 

7 T/ 7 T + 

1.1(1) 

1.2 

1.1 

1.1 

77/7T° 

0.081(13) 

0.099 

0.121 

0.109 

K a s 

0.125(19) 

0.149 

0.172 

0.159 

K" s /hr 

0.123(20) 

0.131 

0.149 

0.138 

A /hr 

0.077(11) 

0.102 

0.076 

0.091 

A IK 

0.63(8) 

0.78 

0.51 

0.66 

K+/K- (NA44) 

1.85(9) 

1.68 

1.78 

1.74 

K+/K* (NA49) 

1.8(1) 

1.68 

1.78 

1.74 

S+/A 

0.188(39) 

0.142 

0.291 

0.191 

(H++H-) 

( A+A ) 

0.13(3) 

0.09 

0.16 

0.12 

Z+/Z- (NA49) 

0.232(33) 

0.223 

0.237 

0.232 

Z+/Z~ (WA97) 

0.247(43) 

0.223 

0.237 

0.232 


0.383(81) 

0.448 

0.522 

0.493 

n~/z~ 

0.219(45) 

0.133 

0.136 

0.135 


Table 4.1: Thermal-model results for the particle ratios measured in Pb+Pb collisions at CERN 
SPS. The experimental data are used in the form compiled in Ref. |E|. The contribution from the 
weak decays is treated in three different ways: i )maximal feeding (third column) means that the 
contributions from all weak decays are included in the multiplicities of the produced hadrons, ii) no 
feeding (fourth column) means that the weak decays have been reconstructed, and the contributions 
from the weak decays are not included, and iii) 50% feeding means that 50% of the contribution 
from the weak decays is included in the final multiplicities. 
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shown in the third column), A a /h a in the second version (no feeding, the fourth column), 
and (A a + A^) / (h~ + ) in the third version (averaged feeding, the fifth column). Other 

ratios are treated in the similar way. In the second column the experimental data are shown 
according to Ref. JT3| . We observe that different treatment of the weak-decay contributions 
leads to slightly different values of T chem , ^f hern , ^ hem and /A em , but the three cases are 
consistent with each other within errors (compare the first two raws of Table [hip . 


The best fit is obtained when the feeding from the weak decays is averaged, which cor¬ 
responds most likely to the experimental situation. In this case: T c h em = 168 ± 3 MeV 
and i^chem = 244 ± 8 MeV. The authors of Ref. |H| obtained T chem = 168 ± 2 MeV and 
Hchem — 266 ± 5 MeV. The agreement of the fitted temperatures is very good, the difference 
of the fitted values of the baryon chemical potential may be caused by a different treatment 
of many branching ratios which are not well known. In the most realistic case (50% feeding) 
our x 2 per one degree of freedom is smaller than one, indicating a good quality of the thermal 
fit, see Fig. |4.1| where the contour plot of y 2 /n is shown. The total number of degrees of 
freedom is 20 in our case. It is smaller by one from the number of ratios included in Ref. |TTJ| , 
since we discard the deuteron measurement. The latter may be disturbed by the inclusion 
of nuclear fragments. 



/Ahem [MeV] 


Figure 4.1: The contour plot of y 2 /n for the ratios measured in Pb+Pb collisions at CERN SPS. 
The numbers at the contours give our values of y 2 /n. Our optimal values of T c h em and are 

denoted by the black circle. The optimal values of Ref. |l3|] are denoted by the triangle. The 
crosses mark the errors. 
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Pb+Pb 

experiment 

quantum 

statistics 

classical 

statistics 

m < 1.8 GeV 

m < rrifi 

Tchem [MeV] 


168+3 

167+3 

171+3 

171+3 

A™ [MeV] 


244+8 

243+8 

247+8 

248+8 

Pchem [MeV] 


62 

62 

62 

61 

Vthem t MeV ] 


-9 

-9 

-9 

-8 

X /n 


0.9 

1.1 

0.9 

1.1 

{Pa ~Pa) IK 

0.228(29) 

0.224 

0.228 

0.228 

0.212 

Pal Pa 

0.055(10) 

0.060 

0.060 

0.061 

0.060 

Pb/Pb 

0.085(8) 

0.081 

0.081 

0.082 

0.083 

A c /A c 

0.131(17) 

0.118 

0.118 

0.118 

0.117 

S-/Ac 

0.110(10) 

0.105 

0.105 

0.115 

0.123 

S+/A c 

0.206(40) 

0.207 

0.207 

0.229 

0.240 

7l~/ tT + 

1.1(1) 

1.1 

1.1 

1.1 

1.1 

r)/i r u 

0.081(13) 

0.109 

0.111 

0.109 

0.114 

iq /*- 

0.125(19) 

0.159 

0.162 

0.156 

0.161 

Ks /h- 

0.123(20) 

0.138 

0.140 

0.136 

0.140 

A /hr 

0.077(11) 

0.091 

0.093 

0.090 

0.092 

A/ATg 

0.63(8) 

0.66 

0.66 

0.67 

0.66 

K+/K~ (NA44) 

1.85(9) 

1.74 

1.73 

1.76 

1.78 

K + /K~ (NA49) 

1,8(1) 

1.74 

1.73 

1.76 

1.78 

S+/A 

0.188(39) 

0.191 

0.192 

0.208 

0.217 

(H++H-) 

0.13(3) 

0.12 

0.12 

0.13 

0.14 

(A+A) 

S+/S" (NA49) 

0.232(33) 

0.232 

0.233 

0.231 

0.226 

S+/S- (WA97) 

0.247(43) 

0.232 

0.233 

0.231 

0.226 

n+/n- 

0.383(81) 

0.493 

0.498 

0.486 

0.471 

Q-/E~ 

0.219(45) 

0.135 

0.134 

0.140 

0.141 


Table 4.2: Thermal-model results for the particle ratios measured in Pb+Pb collisions at CERN 
SPS. The experimental data are taken from Ref. [pRfl . The role of the quantum statistics and the 
mass cut is displayed. 


38 





































4.2 Sensitivity to Classical Statistics and Limited Mass 
Spectrum 


In this Section, we discuss the sensitivity of our results to two different modifications. The 
first modification is simply a replacement of the quantum distribution functions by the classi¬ 
cal distributions (the Bose-Einstein distribution for mesons and the Fermi-Dirac distribution 
for baryons are replaced by the Boltzmann distribution for all hadrons including pions and 
kaons). The second modification is connected with the introduction of a limiting hadron 
mass. In this way we can check how important the tail of the mass distribution of hadrons 
is for the results of our analysis. 


The results obtained with the classical distribution functions and with the limited hadron 

The data are listed in the second column. In the 


mass spectrum are shown in Table [T2 


third column we present our results obtained for exact quantum statistics. They coincide 
with the last column of Table |4.1| , i.e., we average the feeding from the weak decays here. 
In the fourth column we show the results for the classical statistics. The fifth and sixth 
columns show our results obtained with the limited mass spectrum. In the fifth column we 
neglect the feedback from the resonances heavier than 1.8 GeV, whereas in the sixth column 
the maximal mass is equal to the mass of (uiq = 1.672 GeV). In each column we show 
the corresponding values of the optimal thermodynamic parameters and x 2 /n. 


We find that the use of the classical distribution functions leads to very small changes of 
our results. The values of T chem and change only by 1 MeV, and y 2 /n increases from 

0.9 to 1.1. The quantum statistics are not important in the thermal analysis because most of 
particles are very heavy and they can be described well by the classical distributions. Even 
for pions we can use the classical statistics, since they are produced mainly by the decays of 
heavier resonances. The classical features of the hadronic fireball at the freeze-out are also 
reflected in the equation of state, which has a form P = (um + ub)T. Cutting the mass 
spectrum at 1.8 GeV changes the values of T c j iem and fi^ hem by 3 MeV. The cut at rriQ = 
1.672 GeV has a similar impact on T c j iem and Uchemi which remain in agreement (within 
errors) with the exact results. We have checked that the cuts at smaller masses, ~ 1.6 GeV, 
cause already appreciable changes of T c h e m and /x^ em . 


4.3 Thermodynamic Properties of Freeze-Out 

In this Section we present the complete set of the thermodynamic parameters character¬ 
izing the chemical freeze-out in Pb + Pb collisions at CERN SPS. We have again done 
our calculations in three different ways, treating differently contributions from the weak 
decays. As discussed above, in the three considered cases the values of T chem and ^hem 
are consistent with each other. It is interesting to observe, however, that small differences 
in T chem and ^f hem may result in quite large changes of other thermodynamic parameters, 
obtained as the integrals over the distribution functions (|2.2|) . In Table |4.3| we show the 
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Pb+Pb 

Max. feeding 

No feeding 

50 % feeding 

Tchem [MeV] 

164±3 

170±3 

168±3 

Vchem [MeV] 

234±7 

250±8 

244±8 

Rcfcem ^V] 

56 

65 

62 

Rcfcem [MeV] 

-8 

-9 

-9 

Eb [GeV/fm 3 ] 

0.24 

0.35 

0.31 

Em [GeV/fm 3 ] 

0.35 

0.45 

0.41 

e [GeV/fm 3 ] 

0.59 

0.80 

0.72 

P [GeV/fm 3 ] 

0.09 

0.12 

0.11 

p B [1/fm 3 ] 

0.13 

0.19 

0.17 

s [1/fm 3 ] 

4.0 

5.1 

4.7 

tib [1/fm 3 ] 

0.15 

0.22 

0.19 

n M [1/fm 3 ] 

0.40 

0.49 

0.46 

r B [GeV] 

1.6 

1.6 

1.6 

r M [GeV] 

0.9 

0.9 

0.9 

r [GeV] 

1.1 

1.1 

1.1 


Table 4.3: Thermodynamic parameters at the chemical freeze-out in Pb + Pb collisions at CERN 
SPS. Different treatment of the weak decays leads to quite large changes in the estimates of the 
energy density, baryon number density, and other thermodynamic quantities. 


energy density, pressure, baryon number density, entropy, baryon and meson densities, and 
the Cleymans-Rccllich ratios. We find that different treatment of the weak decays causes 
that the thermal-model estimate of the energy density at the chemical freeze-out varies from 
0.6 GeV/fm 3 to 0.8 GeV/fm 3 . Similarly, the baryon density at the freeze-out changes from 
0.13 frn -3 to 0.19 fm -3 . Such changes indicate, that the correct experimental reconstruction 
of the weak decays is necessary in order to have more precise information about the state of 
matter at the freeze-out. The last three raws of Table |4l| show the Cleymans-Rccllich ratio 
calculated for baryons (re), mesons (tm), and all hadrons (r). We observe that r is close to 
unity in the three considered cases, i.e., it is insensitive to the way in which the weak decays 
are treated. We again find that is much larger than r«. 


4.4 Scaling of Masses and Widths at Freeze-Out 


Thermal-model fits show that the temperature at the chemical freeze-out, T c h em , as well as 
the baryon density, are large. As we have just seen in the previous Sections, one typically 
obtains T chem 170 MeV, which is close to the expected critical value for the deconfine- 
ment/hadronization phase transition [fk|. In this situation one may expect that hadron 


properties at the chemical freeze-out are strongly modified by the presence of the hadronic en- 
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vironment. Indeed, such modifications are predicted by different model calculations [28 — 31], 
which helps to explain the low-mass dilepton enhancement observed in the CERES ]53| and 
HELIOS [Q experiments. In this Section we incorporate possible modifications of hadron 
masses and widths into thermal analysis of the particle ratios. We generalize the results of 
Refs. |f43[ , |44| where the problem was studied without refitting thermodynamic parameters. 

Initially, we include only the mass modifications and calculate the particle densities from 
the ideal-gas expression 


rii = 


9i 


p 2 dp 


2n 2 Jo exp [( E i dchem^i dchem^i Rchem^i) /^~ chem ] + e 


where 


E* = \Jp 2 + ( m*Y 


(4.1) 


(4.2) 


is the energy. Of course, in standard thermal-model fits Eq. (|4.1|) is used with the vacuum 
masses, m* = The in-medium masses, m*, may depend on temperature and density in a 

complicated way. In order to explore possible different behavior of in-medium masses and, 
at the same time, keep simplicity, we do our calculations with the meson and baryon masses 
rescaled by the two universal parameters, Xm and Xg, namely 


m M ~ x m ~ x b bib- 


(4.3) 


An exception from this rule are the masses of pseudo-Goldstone bosons (n, K, rf) which we 
keep constant. This is in agreement with explicit model calculations incorporating chiral 
symmetry, e.g. |55], We note that the use of Eq. (0) is valid when the in-medium 
hadrons can be regarded as good quasi-particles. A thermodynamically consistent approach 
has been constructed so far only for the lowest multiplets of hadrons 0- At SPS energies, 
however, it is crucial to include all hadrons with masses up to (at least) 1.6 GeV. For 
such a complicated system, a thermodynamically consistent approach is not available. As 
in the standard approach, Eq. ([44]) is used to calculate the primordial density of stable 
hadrons and resonances at the chemical freeze-out. The final (observed) multiplicities receive 
contributions from the primordial stable hadrons, and from the secondary hadrons produced 
by decays of resonances after the freeze-out. Ref. [[40| is used to determine the branching 


ratios, which we keep unchanged. We neglect the finite-size and excluded volume corrections 
which do not affect the particle ratios. 

For fixed values of x M and Xb we fit the temperature, T chem , and the baryonic chemical 
potential, pfhemi by the minimization the expression 


X 2 (xm,xb) = X ) 


[ R 


exp 

k 


td therm 

-*4 


(%i X B) 


-i 2 


k =1 


crt 


(4.4) 


where Rj xp is the kth measured ratio, cr*, is the corresponding error, and RjJ ierm (x m, x b) is 
the same ratio as determined from the thermal model with the modified masses. Similarly to 
the standard case with unchanged masses, m* = mi, the strange chemical potential, /nf hem , 
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and the isospin chemical potential, Hc hem , are determined from the requirements that the 
initial strangeness of the system is zero, and that the ratio of the baryon number to the 
electric charge is the same as in the colliding nuclei, see Eqs. (|2.5|) and (|2.6|). 



0.8 0.9 1 1.1 1.2 
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Figure 4.2: Dependence of y 2 , and the fitted values of the temperature and the baryon density on 
the scale parameter x. Plot represents Pb + Pb collisions at SPS energies. Solid lines: all hadron 
masses (except for Goldstone bosons) are scaled with xm = xb = x. Dashed line: only meson 
masses are scaled, xm = x,xb = 1. Dotted lines: only baryon masses are scaled, xb = x,xm = 1 . 
The nuclear saturation density po = 0.17 fm -3 . 


In Fig. |4.2| we plot our results obtained for the experimental ratios measured in Pb + 
Pb collisions at SPS. In the case xm — %b — 1 we recover our results discussed in detail in 


Sections |4.1| - |4.3| . In Fig. [4. 2| a) we give our values of y 2 . One can observe that a small 
decrease of the meson and baryon masses, % = xb ~ 0.9, leads to a slightly better fit with 
the corresponding smaller values of the temperature and the baryon density, as shown in 
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Figs. [h^ b) and fO| c). It would be premature, however, to conclude that the masses drop. 
The values of y 2 for the solid line are increased by 25% compared to the minimum in the 
range 0.75 < x < 1.05, which clearly is the allowed range. We thus conclude that moderate 
dropping of hadron masses, say by 20%, does not spoil the quality of thermal fits. On the 
contrary, larger dropping or growing of the masses result in a significant increase of the y 2 
values. 


With the modified masses the thermodynamic parameters characterizing the fits change. 
For example, if we rescale both meson and baryon masses (except for Goldstone bosons) in 
the same way, x = xm = %b, the temperature and the chemical potentials are to a very 
good approximation also rescaled by x. This follows from the fact that we study a system of 
equations which is invariant under rescaling of all quantities with the dimension of energy. If 
we allowed also for the changes of the masses of the Goldstone bosons, the thermodynamic 
parameters would scale exactly as T chem (x ) = x T chem (x = 0) and p c hem(x ) = x p c hem{x = 0). 
In this case y 2 remains constant, independently of x. For fixed values of the Goldstone-boson 
masses, the scale invariance is broken, y 2 varies with x, as shown in Fig. |4.2| a), and the 
results are non-trivial. 


To account for finite in-medium widths, T*, of the resonances we generalize Eq. O to 
the formula [57-60] 


m = / dM 2 


Ml 


CXJ 

I dp 


m* r* 


x 


nN (M 2 - (m *) 2 ) 2 + 


p 


2vr 2 exp [(V M 2 + p 2 - ^ hern B, t - n s chem Si - pi hern Ii) /T chem 


± 1 


where N is the normalization of the relativistic Breit-Wigner function, 


Ar 1 1 

— —|— arctan 

2 7T 


(K) - Mji/m 


i. 


(4.5) 


(4.6) 


The integral over M 2 is taken to start at the threshold Mq corresponding to the dominant 
decay channel. In the limit T* —> 0 Eq. ( |4.5D obviously reduces to formula (|4.1|) . 

In order to analyze the effect of broadening of hadron widths we introduce the parameter 
y in such a way that 

r; = y Ti. (4.7) 

Here T* are the vacuum widths, hence the case y — 1 corresponds to the physical widths as 
measured in the vacuum, and the case y = 0 represents the situation when the widths are 

we show the results of our 


neglected (our previous analysis based on Eq. ( |4.1|) ). In Fig. 
fitting procedure. We observe that the inclusion of the vacuum widths does not change the 
value of y 2 , and the values of T chem and p^ em . An increase of the widths by a factor of 2 
has also little effect. Only for larger modifications of the widths the fit gets worse. 

In conclusion we state that the thermal analysis of particle ratios, measured in Pb + Pb 
collisions at CER.N SPS, allows for moderate dropping of hadron masses (~ 20%). This does 
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Figure 4.3: Dependence of y 2 , and the optimal values of the temperature and the baryon density 
on the scale parameter for the widths, y. 


not spoil the fits, which remain of similar quality as those obtained without modifications. 
Larger dropping of the hadron masses or growing of the masses are not likely. Scaling of 
hadron masses results in modifications of the thermodynamic parameters for which the fits 
are optimal. In particular, lowering of all the masses leads to a smaller values of T c h em and 
Pchem- This might be a desired effect, since T c h em ~170 MeV is large and may correspond 
to quark-gluon plasma rather than to a hadron gas. Our study of the modifications of the 
hadron widths shows that they have small impact on the ratios. 
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Chapter 5 

Au + Au Collisions at BNL RHIC 


The heavy-ion program at CERN SPS delivered several very interesting results, but there was 
no clear discovery of the new physical phenomena, such as the desired deconfinement phase 
transition. The newest heavy-ion machine, the Relativistic Heavy-Ion Collider constructed 
at BNL, was designed to accelerate heavy ions at energy a/s = 200 GeV. In the year 2000, 
during the first run of RHIC, the maximal energy of /s = 130 GeV was achieved. This 
energy exceeds the CERN SPS energy by one order of magnitude [J so the new phenomena 
were expected to occur. The first new RHIC data indicate indeed at several new features 
of the collision process such as: much higher particle multiplicities, increased production of 
antiparticles, and lower baryon number in the central region. Nevertheless, the overall picture 
of the collision follows the pattern established from the studies of heavy-ion collisions at lower 
energies. Further systematic analysis of the just incoming data is necessary to extract more 
information. 


5.1 Thermal Analysis of Particle Ratios 

Table |5d] presents our fit to the particle ratios measured at RHIC during its first run at a/s = 
130 GeV. We stress that exactly the same version of the thermal model has been used in this 
fit, as that used in our previous studies of Si + Au, S + Au, and Pb + Pb reactions. In our 
present calculation, the identical ratios measured by different groups are treated separately 
in the definition of y 1 2 (number of points n = 16). In this way the measurements done by 
different groups enter independently, and converging experimental data have a larger weight 
in the statistical analysis. Very similar results are obtained, however, if we first average the 
results of different groups to obtain the most likely value for each considered ratio. This fact 
shows consistency of the experimental measurements done by different groups, see the third 
column in Table |5.1| . 

1 The energy of 158 GeV per nucleon in the lab corresponds to /s = 17 GeV per nucleon pair in the 

center-of-mass system. 
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Au+Au 

TM 

experiment 

Tchem [MeV] 

165±7 


tchern [MeV] 

41±5 


Pchem [M e V] 

9 


^chem [M e V] 

-1 


X/ n 

0.97 


n~/ ti + 

1.02 

1.00 ±0.02 

61 

, 0.99 ±0.02 

52 


p/tt~ 

0.09 

0.08 ±0.01 

M 


K~/K+ 

0.92 

0.88 ±0.05 

0.91 ±0.09 

M 

|6l 

, 0.78 ±0.12 

, 0.92 ±0.06 

|65 

|62 


K~f 7 T “ 

0.16 

0.15 ±0.02 

64 


K'Jh- 

0.046 

0.060 ±0.012 

HM 


K'Jh~ 

0.041 

0.058 ±0.012 

34|,|66| 


P/P 

0.65 

0.61 ±0.07 

0.60 ±0.07 

61 

, 0.54 ±0.08 
, 0.61 ±0.06 

|65 

52 


A/A 

0.69 

0.73 ±0.03 

H 


5/5 

0.76 

0.82 ±0.08 

M 



Tabic 5.1: Thermal fit of the particle ratios measured at RHIC at y/s = 130 GeV. 


Our optimal value of T c h ern obtained from the analysis of the particle ratios equals 165 ± 7 
MeV. It is a very interesting fact that T c h em for RHIC agrees well with T c h em = 168 ± 3 MeV 
found in our fit for Pb + Pb collisions at CERN SPS. The growth of the energy of the 
colliding nuclei by one order of magnitude does not lead to creation of a hotter system, only 
the baryon chemical potential is significantly lower than that found at CERN SPS. Clearly, 
the collisions at RHIC energies are more transparent than the collisions at CERN energies. 
We have also calculated other characteristics of the freeze-out. In particular, we find the 
energy density e = 0.5 GeV/frn 3 , the pressure P = 0.08 GeV/fm 3 , and the baryon density 
Pb = 0.02 fnu 3 . Here again, one can observe that the energy density is not higher than that 
found at CERN SPS, only the net baryon number tends rapidly to zero. Our calculation 
confirms the Cleymans-Redlich conjecture |14| that the energy per hadron at the chemical 
freeze-out is 1 GeV (our approach yields almost exactly r = 1.0 GeV). We observe, however, 
that the average baryon energy is much larger than the average meson energy: r B = 1.6 GeV 
and tm = 0.9 GeV. We emphasize that exactly the same r B and have been extracted 
from Pb + Pb collisions at CERN SPS, see Table [Oj Also the values of r B and r M extracted 


from the sulphur collisions, see Table |3.6| , agree well with these two values. Further work 
is needed to explain such universal behavior. In addition, in our thermal approach we find 
that the ratios A/A and 5/5 are practically unaffected by the weak decays, since the latter 
contribute in the same way to the abundances of baryons and antibaryons (note the very 
small value of our Pchem)- 
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Figure 5.1: The contour plot of our y 2 /n. Our result (black circle), the result of Ref. [|E]] (black 
square), and the fit of Ref. |2(J (black triangle) are all shown with the corresponding errors. The 
numbers at the contours indicate the values of x 2 / n - 


In Fig. [1 .1 1 we show y 2 /n as a function of the temperature, T chem , and the baryon chemical 
potential, R^ em . An interesting feature of the plot is that it shows a characteristic valley of 
the optimal values of the thermodynamic parameters. The shape of this valley indicates that 
the quality of the fit does not change much if we moderately increase or decrease both T chern 
and Uchem- We note that a similar valley can be seen in Fig. |47T| , where the values of x 2 / n 
were plotted for the case of Pb + Pb collisions at CERN SPS. While our work was nearing 
completion [^1], a fit by Braun-Munzinger, Magestro, Redlich and Stachel was presented: 
Tchem = 174 ± 7 MeV and /i^ em = 46 ± 5 MeV PD| . We note that our T c h em is 9 MeV lower 

|, and 25 MeV lower than 190 MeV obtained in a similar fit by Xu and 
Still, the results of the three calculations are consistent within errors. This fact 


than T c hem of Ref. 
Kaneta |19[ 


is displayed in Fig. p~T| , where our result (black circle), the fit of Ref. ]2(J (black triangle), 
and the fit of Ref. |19| (black square) are shown all together with the corresponding errors. 
Especially, the results of Refs. 


and Ell] are close and line up along the valley of the 


optimal parameters. 
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Figure 5.2: Dependence of y 2 on the scale parameter x. In this case xb = xm = x, he., the masses 
of baryons and mesons are rescaled in the same way (except for the pseudo-Goldstone bosons). 


Similarly to the case of Pb + Pb collisions at CERN SPS, we have studied the influence 
of possible in-medium mass modifications on the thermal analysis of the ratios measured at 
RHIC. In Fig. |5.2| we show the dependence of y 2 on the scale parameter x (in this case the 
mass modifications of baryons and mesons are assumed to be equal). We observe that the 
values of y 2 are slightly lowered for x < 1. This is analogous behavior to that observed in Pb 
+ Pb collisions studied in Sect. 4.4. Again, the effect is not dramatic (y 2 decreases by 20% if 
x changes from 1 to 0.8), so the indication for dropping masses is rather weak. On the other 
hand, we see that the changes of the masses can be included in the thermal analysis, yielding 
the fits of the same quality as that found in the calculations without any modifications. In 
the present case, however, we see that larger dropping of the masses is allowed, as compared 
to the CERN-SPS case discussed before. We have also analyzed Si + Au collisions at AGS, 
and S + Au collisions at SPS. In these two cases y 2 has a flat minimum at xm ~ xb ~ 1. 


5.2 Pressure of Hadron Gas at Freeze-Out 

In Fig. |5.3| we plot the pressure of a hadron gas as a function of the temperature and the 
baryon chemical potential. We show the neighborhood of the optimal parameters: T c h em = 
165 MeV and /i^ em = 41 MeV. The exact result shown in Fig. [54], obtained as a sum of all 
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partial pressures given by Eq. ( |2.29|) , can be very well approximated by the formula 


Mmax 

Phg(T,h b ) = — f p M (m)m 2 T 2 K 2 (^j dm 


o 

Mmax 


^ cosh ^ ^ 


B' 


m 2 T 2 Jl 2 


7T 


T 




!27T 2 


h 


S' 


m 


pM(m) + 2 cosh I — I p B (m) 


dm 


7Ti 2 T 2 K 2 dm. (5.1) 


In Eq. (|5.1|) the sum over all hadronic states is replaced by the convolution of the classical 
formula for the partial pressure, Eq. (|2.37|) , with the mass spectrum for baryons and mesons, 
Eqs. (|2.22| ) and ( |2.23| ). The maximal mass included in the integral is taken to be the same 
as that used in the fits of the spectrum |47|], M max = 1.8 GeV. We note that the strange 


chemical potential as well as the isospin chemical potential have been neglected in Eq. ( |5.1| ) . 
Our fit shows that these two potential are very small at RHIC. For T = 165 MeV and 


p B = 41 MeV, the pressure calculated from Eq. (|5.1|) agrees very well with the exact result 


P = 0.08 GeV/fm 3 . The very weak dependence of the pressure on the baryon chemical 
potential, as displayed in Fig. [5.3| . follows also from Eq. (|5TT1) , since in the considered region 
one Ends cosh (/i s /T) ~ 1. 



Figure 5.3: Pressure of a hadron gas which includes all resonances consisting of up, down, and 
strange quarks. 
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Figure 5.4: Our optimal thermodynamic parameters. The fitted values of the temperature for Pb 
+ Pb collisions at CERN SPS and for Au + Au collisions at BNL RHIC are consistent with each 
other, and close to the critical temperature predicted by the QCD simulations on a lattice [49]. The 
lattice results (grey bands) correspond to two different versions of the calculations (the number of 
massless flavors, Nf, is taken to be either 2 or 3). The vertical size of the bands denotes the error 
of the lattice calculations. We show also our optimal parameters describing the collisions of lighter 
nuclei. Note that in this case the quality of the fits is worse, as discussed in detail in Chapter 3. 


It is constructive to compare the pressure of the hadron gas at freeze-out to the pressure 
of the ideal gas of quarks and gluons at the same temperature and baryon chemical potential. 
To do so we use the formula 


Pqgp(T , fx q ) — - 


7T~ 


16 30 + 6Nf 


m Ti+ pF 2 + g? 


— B, 


(5.2) 


where Nf is the number of quark flavors and B is the bag constant. For T = 165 MeV, 
H q = fi B /3 = 41/3 MeV (the quark chemical potential is one third of the baryon chemical 
potential), B l P = 200 MeV, and Nf = 3 one finds P = 0.3 GeV/frn 3 . This result shows 
that the pressure of the ideal gas of quarks and glnons is much larger than the pressure of 
a hadron gas (we recall that onr calculation gives P = 0.08 GeV/frn 3 ). If the equation of 
state (|5.2|) was realistic, we would deal with the plasma rather than with a hadron gas at the 
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temperatures as high as T — 165 MeV. More realistic equations of state of the plasma include 
interactions of quarks and gluons which lower the pressure of the plasma and increase the 
temperature of the possible phase transition. The most reliable predictions concerning the 
phase transition are obtained by the numerical studies of QCD on a discretized space-time 
lattice j49|. The most recent values of the critical temperature are: T c = 173 ± 8 MeV for 


Nf = 2, and T c = 154±8 MeV for Nf = 3. It is intriguing that the thermal fits of the particle 
ratios yield the temperature, which is very close to the critical temperature inferred from 
the lattice simulations. The closeness of these two temperatures suggests that the chemical 
content of the hadronic fireball is established just after hadronization phase transition or, 
some authors argue, during the hadronization process itself. In the latter case the thermal 
distributions of hadrons have little in common with rescattering. Note, however, that our 
fit is consistent (within errors) with a temperature smaller than T c . 



Figure 5.5: Baryon density of a hadron gas in the vicinity of the freeze-out point: T c h em = 165 
MeV and p^ hem = 41 MeV. 


We close this Section with the two additional plots. In Fig. |5.5| we show the baryon 
density of a hadron gas in the region close to our optimal thermodynamic parameters. In 
the considered region the baryon density is very low. It is much smaller than the baryon 
saturation density p 0 — 0.17 frn -3 (we note that for Pb + Pb collisions at CERN SPS the 
fitted baryon density was close to p 0 ). Finally, in Fig. ^6] we plot the Cleymans-Rccllich 
ratio. One can observe that r is close to unity (0.8 < r < 1.2) in the rather vast range of 
the temperature and the baryon chemical potential. 
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Figure 5.6: Cleymans-Redlich ratio for a hadron gas. 
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Chapter 6 
Conclusions 


In the end of the paper we summarize the main findings of our work: 


• Thermal description of the particle ratios, based on the assumption of the full chemical 
and thermal equilibrium , works well for Pb + Pb collisions at CERN SPS and An + An 
collisions at BNL RHIC. At smaller beam energies and for smaller sizes of the colliding 
systems this version of the thermal model reproduces only a qualitative behavior of 
the data. We stress that this observation follows from the study which uses exactly 
the same framework for all reactions. 

• The optimal values of the temperature inferred from the thermal analysis of the ratios 
in Pb + Pb collisions at CERN SPS and An + Au collisions at BNL RHIC agree. The 
more energetic collisions at RHIC do not produce a hotter hadronic system. The fitted 
value of the temperature ( T c h em ~ 165 MeV for RHIC and T c h e m ~ 168 MeV for SPS), 
is very close to the critical temperature obtained from the lattice simulations of QCD. 

• In the case of Pb + Pb collisions at CERN SPS energies, the correct reconstruction 
of the weak decays is essential for the estimates of the energy density and the baryon 
number density characterizing the freeze-out. On the other hand, at RHIC the role of 
the weak decays is less important. 

• Our calculations confirm the Cleymans-Redlich observation that thermal conditions at 
freeze-out, for different colliding systems at different energies, correspond to an average 
hadron energy 1 GeV. Moreover, we have found that averaged baryon energy is 1.6 
GeV and the average meson energy is 0.9 GeV for both Pb + Pb collisions at CERN 
SPS and Au + Au collisions at BNL RHIC. This new observation is parallel to the 
argument that two different Hagedorn temperatures are required to describe the mass 
spectra of baryons and mesons. 
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• Possible modifications of the hadron masses and widths were incorporated into a ther¬ 
mal analysis of the particle ratios. In the case of Pb + Pb collisions at CERN SPS, 
we have found that moderate, up to 20%, dropping of the masses does not spoil the 
quality of the fits. Larger dropping or growing of the masses are not likely, since they 
result in a significant increase of the y 2 values. We have also showed that the increase 
of hadron widths by less than a factor of two does not affect the thermal-model fits. 
Similar behavior is found for Au + Au collisions at BNL RHIC. In this case, however, 
larger dropping of the masses is allowed. 

• We have checked, for a variety of colliding systems, that the quantum statistics (Bose- 
Einstein and Fermi-Dirac) can be very well approximated by the classical Boltzmann 
statistics. This means that the hadron gas at the freeze-out behaves like a classical gas 
with the equation of state P = nT. In addition, the use of the Boltzmann distribution 
function implies that the excluded-volume corrections do not change the values of the 
fitted optimal values of T chem and fifhem (provided the eigenvolumes of baryons and 
mesons are equal). 

• Finally, we emphasize that our fits may be treated as a first step in more complex 
investigations of hadron production in ultra-relativistic heavy-ion collisions. A com¬ 
bination of the thermal approach with the hydrodynamic expansion can be used to 
study other observables such as: transverse-momentum spectra, rapidity distributions, 
elliptic flow, or HBT correlation radii. 
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